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1. Introduction 



Polarized deep-inelastic scattering (DIS) has attracted considerable theoretical Q and 
experimental interest since the EMC experimental results in 1988 Q showed that the 
proton spin structure is subtler than naive parton expectations might suggest. A first 
generation of experiments with electron and muon beams, and proton, deuterium and 
neutron (i.e. 3 He) targets has provided us with information on individual polarized parton 
distributions which give important clues on the relevant underlying theoretical issues. Sev- 
eral crucial pieces of information, however, are still missing. Specifically, current inclusive 
DIS experiments only measure the C-even combination Aq + Aq of quark distributions. 
Moreover, they are only very weakly sensitive to the size of the strange quark distribution. 
Therefore, a full flavor separation is very difficult, and a separate measurement of the quark 
and antiquark components is impossible. Furthermore, there are indications MM that the 
polarized gluon distribution is large, but uncertainties are still sizable. Several of these 
pieces of information are crucial not only in reconstructing the full spin structure of the 
nucleon, but also (perhaps more interestingly) in deciding among the various theoretical 
scenarios which have been proposed to understand the somewhat puzzling picture which 
the original EMC results and their subsequent refinements point at. In particular, it is 
crucial to determine accurately the size of the strange contribution || and of the polarized 
gluon 0. 

In the unpolarized case, the most precise available information on the strange 
quark distribution and the quark-antiquark separation comes from charged-current deep- 
inelastic scattering ||, which, being mediated by a charged current, gives access to the 
flavor structure of the target. In the polarized case, the luminosity of existing neutrino 
beams requires targets whose typical sizes are of the orders of several meters, and therefore 
cannot be polarized. In principle, charged-current events can be a significant part of the 
cross section even with electron or muon beams if the energy is high enough. Indeed, 
charged-current unpolarized cross sections and structure functions have been measured at 
HERA U, and the corresponding polarized asymmetries have been studied fl(| as a way 
of gaining information on polarized parton distributions if the polarized option were to be 
available at HERA. 

It is several years now since muon storage rings have been proposed as a powerful 
new way to reach high energies and luminosities in a colliding-beam facility [^TJ,|TJ[. It 
has since been recognized that exciting physics can also be obtained with the highly fo- 
cused neutrino beams arising from the decays of muons along straight sections of the 



accumulator |T3[ . In addition to the obvious applications for studies of neutrino properties 



(masses and mixings) using far-away detectors (for a recent detailed study see Ref. Hl4fl ) 



DIS experiments operating close-by downstream the muon ring could provide significant 
contributions to several topics |15],[16[] : from the study of the nucleon structure to accurate 
measurement of the parameters of weak and strong interactions, from the determination of 
CKM matrix elements to studies of heavy quark decays, and more. Reference []15[ pointed 
out the potential for measurements of unparalleled precision of both unpolarized and po- 
larized neutrino structure functions, leading to an accurate decomposition of the partonic 
content of the nucleon in terms of individual (possibly spin-dependent) flavor densities. A 
more quantitative investigation of polarized DIS with neutrino beams is therefore called 
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for, both in order to assess the physics potential of these experiments, and also to decide 
among other possible options, such as polarized HERA. A first study in this direction has 
recently appeared in JT7J. 

Several studies of polarized neutrino DIS have been presented in the literature. In 
particular, the parton model formalism has been presented in Ref. |18|], while its relation 



to perturbative QCD has been worked out through the operator-product expansion at 



leading order in Refs. [19-21]. The next-to-leading order coefficient functions were com- 



puted several years ago [22 , and have been recently generalized to the case where heavy 



quark masses are retained |23[ ; next-to-leading order polarized anomalous dimensions were 
determined for the full set of operators relevant to both neutral- and charged-current DIS 
recently |[24| . Even though all the theoretical tools which are needed in order to perform 



a full next-to-leading order analysis are thus essentially available in the literature, no such 
study has been performed yet. 

In this paper, we wish to give a first systematic next-to-leading order treatment of 
charged-current polarized DIS. The purpose of this work is to provide a set of common 
theoretical tools, as well as some benchmark estimates for future work. We will start 
(Sect. 2) by recalling the structure of the polarized asymmetries in charged-current DIS 
in terms of structure functions, review their leading-order expressions in terms of parton 
distributions, and the scale dependence of the parton distributions. We will then recall 
how charged-current lepton and antilepton DIS on proton and neutron targets allows a full 
flavor decomposition of the nucleon spin content. In Sect. 3, we will discuss the full next- 
to-leading order formalism for structure functions and parton distributions (this section 
can be skipped by those who are only interested in the phenomenological analysis). The 
results presented here have been incorporated in an upgrade of the ABFR || next-to- 
leading order evolution code for polarized parton distributions, by including the option of 
evolving charged-current structure functions and fitting them to experimental data. We 
will next (Sect. 4) discuss the constraints that the elementary requirement of positivity 
of cross-sections imposes on the otherwise unknown C-odd combinations of quark and 
antiquark distributions Aqi — Aqi, and how they limit the set of possible expected results 
for these quantities. Then, we will examine the current theoretical expectations for the 
strange and quark-antiquark content of the nucleon within several theoretical scenarios 
which have been suggested in order to understand current information on the nucleon 
spin. 

We will turn to phenomenology in Sect 5, where we will discuss the precision on the 
determination of structure functions which is expected at a future muon storage ring fa- 
cility: in particular, we will review the expected characteristics of the neutrino beam, and 
present detailed estimates of the expected percentage error on structure function determi- 
nations as a function of the kinematic variables x and Q 2 . We will then (Sect. 6) study 
how such a determination of polarized structure functions would improve our knowledge of 
parton distributions: we will generate pseudo-data according to the theoretical scenarios 
of Sect. 4, distributed according to to the error estimates of Sect. 5, and use the evolution 
code discussed in Sect. 3 to perform next-to-leading order fits of structure functions to 
these data. We will thus be able to assess the error on individual parton distributions 
which would be obtained thanks to this kind of experiment. 
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2. Flavor content and scale dependence of polarized structure functions 



Polarized structure functions for charged-current DIS have been computed in the 
parton model in Ref. \TB\, and at leading order in perturbative QCD with operator methods 
in Refs. [p0| , |21[| . Here for completeness we briefly review the derivation of the leading order 
expression of these structure functions from the QCD improved parton model, the evolution 
equations which govern their scale dependence, and then review how it is possible to use 
this information to disentangle the polarized flavor and antiflavor content of the target. 

2.1. Structure functions and parton distributions 

We define structure functions in terms of the hadronic tensor W^ u as follows: 
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where q^, p^ and s M are respectively the momentum of the incoming virtual gauge bosons, 
and the momentum and spin of the incoming nucleon, and 
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The proton spin vector is normalized as s 2 = 
hadronic tensor is defined as |25] 

w ^ = hj d4zei9Z {p,sl ^ w ' J - (0) ] lp > s) - 



2 , where m is the hadron mass. The 



(2-3) 



Notice that the definition with <-> v is sometimes (e.g. in []18[]) adopted; this corresponds 
to changing the sign of the antisymmetric part of the tensor. The definitions of the un- 
polarized structure functions Fj and of the polarized parity-conserving structure functions 
<7i, g 2 are essentially standard, while a variety of different conventions have been adopted 
in the literature for the polarized parity-violating structure functions #3, #4 and g$. Here 



we adopt the same convention as Ref. pi 



We now specialize to the case of charged-current scattering. The total cross section 
for charged-current DIS on unpolarized targets is given by 



d 1 c x t(x 1 y,Q< 
dxdy 
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where A^ = ±1 is the helicity of the incoming lepton (—1 for a neutrino or electron, +1 for 
an antineutrino or positron). The explicit expressions of the various structure functions in 
eq. ( |2.4p depend only on the charge of the exchanged W. and on the target. The polarized 
cross-section difference 

Act = er(A p = —1) — er(A p = +1), (2.5) 
where A p = ±1 is the proton helicity, is given by 
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The results in eqs. (|2.4| ) and ( |2.6| ) agree with those given in Ref. 
factor of two larger than the results of Ref. 
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J, while they are by a 
The definition eq. ( |2.5|) corresponds to 
the difference of antiparallel minus parallel spins of the incoming particles for an incoming 
lepton, but parallel minus antiparallel spins for an incoming antilepton. 

Equation ( |2.6| ) shows that for longitudinal nucleon polarization the contributions of the 

2 

structure functions g 2 and g 3 to the cross-section are suppressed by powers of These 
structure functions give an unsuppressed contribution to the cross-section for transverse 

2 

polarization, but in such case the polarized cross section difference itself vanishes as ^ ~~ * 
0. Henceforth, we will systematically neglect all power suppressed contributions; therefore, 
we will not discuss these structure functions further. It should be pointed out, however, 
that at the neutrino factory it will be possible in principle to extract the power-suppressed 
contributions directly from the data, by studying the Q 2 dependence. 

Because the same tensor structures appear in the spin-dependent and spin- 
independent parts of the hadronic tensor eq. ( |2.1|) in the |p — > limit, the cross section 
difference eq. ( |2.6|) is obtained from the total cross-section eq. ( |2.4|) replacing 



-05, 



-94, 



2<7i 
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and multiplying by a factor two due to the fact that the total cross-section is an average 
over initial state polarizations. In particular, it is easy to derive a polarized analogue |26 



of the Callan-Gross relation [27], which follows from the observation that the quark-gluon 
coupling conserves helicity when all masses are neglected, so the hadronic tensor eq. ( |2.1|) , 
if computed at leading order, must vanish when contracted with a longitudinal polarization 
vector e^. Since can be written as a linear combination of p and this condition implies 
PtiPvW^ 1 ' = 0, which using the expression eq. ( |2.1|) becomes 
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Therefore at leading order 

g 4 (x 1 Q 2 )=2xg 5 (x,Q 2 ) 1 (2.9) 

and, at leading twist, there are only two independent polarized structure functions, namely 
gi (parity conserving) and gs (parity violating), respectively analogous to the unpolarized 
structure functions F± and F3. We will see in Sect. 3.1 that also beyond leading order there 
are only two independent structure functions, despite the fact that eq. (|2.9|) is violated. 

The leading-order expression of g± and g$ are straightforwardly found by computing 
the hadronic tensor for scattering off a free quark ||18|| : 



g? (x, Q z ) = Au(x, Q z ) + Ad(x, Q z ) + Ac(x, Q z ) + As(x, Q< 
gf~ (x, Q 2 ) = Au(x, Q 2 ) + Ad(x, Q 2 ) + Ac(x, Q 2 ) + As(x, Q 2 ) 
gf + (x, Q 2 ) = Au(x, Q 2 ) - Ad(x, Q 2 ) + Ac(x, Q 2 ) - As(x, Q 2 ) 
gf~ (x, Q 2 ) = -Au(x, Q 2 ) + Ad(x, Q 2 ) - Ac(x, Q 2 ) + As(x, Q 2 ) 



(2.10) 



Below charm threshold only the Cabibbo-suppressed part of the As contribution survives, 
because the Cabibbo-enhanced transition would require production of a c quark in the 
final state. A Ac contribution is in principle possible even below charm threshold in 
case intrinsic charm |28[ is present, since it only requires production of a strange quark. 
Similarly, the Cabibbo-suppressed part of the down contribution vanishes. Hence, below 
charm threshold the leading-order expression of the structure functions is 

gf + (x, Q 2 ) = Au(x, Q 2 ) + cos 2 6 c Ad(x, Q 2 ) + Ac(x, Q 2 ) + sin 2 6 c As(x, Q 2 ) 

gY~ (x, Q 2 ) = Au(x, Q 2 ) + cos 2 8 C A d(x, Q 2 ) + Ac(x, Q 2 ) + sin 2 8 c As(x, Q 2 ) 

, (2.11) 
gf (x, Q 2 ) = Au(x, Q 2 ) - cos 2 9 c Ad(x, Q 2 ) + Ac(x, Q 2 ) - sin 2 6 c As{x, Q 2 ) 

gf~ (x, Q 2 ) = -Au(x, Q 2 ) + cos 2 6 c Ad(x, Q 2 ) - Ac(x, Q 2 ) + sin 2 6 c As(x, Q 2 ), 

where 9 C is the Cabibbo angle, and Ac is the intrinsic charm distribution. The remaining 
mixing angles in the CKM matrix can be taken to vanish for all practical purposes, so a 
contribution from Ab will only be present above the top threshold; it will then have the 
form which straightforwardly follows from repeating the pattern of previous generations. 
Equation ( |2.1U| ) is to be contrasted with the familiar case of virtual photon scattering, 
where 



2 — - J (2.12) 



g2 (x,Q 2 )=0. 

2.2. Evolution equations 

The scale dependence of parton distributions is determined by the Altarelli-Parisi 
evolution equations [p9| , |30[] . Define the C-even and C-odd quark distributions 

Aqf = A qi ±Aq u (2.13) 
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the singlet combination 

n f 



AS ± = ^A 9 f; AE = AE + (2.14) 



i=i 



and the nonsinglet combination 



Ag£ s± = (Ag±-Ag±), (2.15) 
The evolution equations have the form 
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eft 2tt s ' 
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where t = log (Q 2 /A 2 ), <g> denotes the usual convolution with respect to x, and all splitting 
functions admit a perturbative expansion of the form AP(x, a s ) = AP^ -* + ^ APW + . . .. 

Notice that the gluon distribution cannot mix with the C-odd distributions because 
charge conjugation implies that 

Pqg = Pqg'i Pgq = Pgq- (2-17) 

Furthermore, because all flavors are equivalent when masses are neglected, it is only neces- 
sary to distinguish flavor-conserving (diagonal) splittings P qq = P qiqi and flavor-changing 
(non-diagonal) splittings P qq D = Pq iqj with i ^ j. Finally, also by charge conjugation, 

Pqq = Pqqi Pqq = Pqq'i A "^W = ^Pqq'i A "^W = ^Pqq- (2-18) 

It then follows immediately that, in general 

AP± S = (AP£ - AP™) ± (AP* - AP™) 
AP± = (AP q D q + (n f - l)AP q N q D ) ± (AP q D q + (n f - l)AP q N q D ) 



(2.19) 



similarly to the unpolarized case (with P — > AP). 

Specializing now to the leading-order case, note that, because the quark-gluon in- 
teraction conserves helicity, flavor, and baryon number, at leading order there is only one 
quark-quark splitting function 

^2= A ^£=^°> (2-20) 



while 

t 

<?<? ii 



(2.21) 

pND, LO _ a pND, hO _ n v ' 

qq QQ 



It follows that at leading order all quark (polarized and unpolarized) distributions evolve 
with the same splitting function: 



AP~ = AP-> S L ° = AP+> s LO = AP^ LO = (2.22) 

2.3. Flavor decomposition 

Assuming the availability of neutrino and antineutrino beams, and the capability of 
measuring independently g\ and g$, we have four different linear combinations of individual 
polarized parton distributions from charged-current scattering, on top of the usual one (or 
two, if proton and neutron targets are available) from neutral-current scattering. Using 



the leading order expressions of the structure functions, eq. (|2.10| ), we get 

\ ( 9 r - 9 r) Ac; \ + 9 r) a c; ^ 

\ ( 9 V + - 9f ) = Arf + As- \ [gf- + gf~) = Ad+ As. (2.24) 



Below charm threshold, in the absence of intrinsic charm, and neglecting Cabibbo- 
suppressed contributions, the Ac and As contributions decouple, and all four light flavors 
and antiflavors can be extracted directly from eqs. ( |2.23 ), (|2.24 ). Above charm threshold, 



or even below threshold if intrinsic charm is present, a direct measurement of all indi- 
vidual parton distributions is only possible if further independent linear combinations are 
experimentally available. 

For a nucleon target, further independent linear combinations of parton distributions 
can be determined by using both proton and neutron (or deuterium) targets: by isospin, 
all structure functions for a neutron target are expressed in terms of parton distributions of 
the proton by interchanging u and d in the expressions for a proton target. This then gives 
us four new combinations, obtained interchanging u and d in eqs. (|2.23| ), (|2.24|) . Combining 
these with the combinations eqs. ( |2.23 ),( 2.24|) , however, only gives us six independent linear 



combinations. A convenient choice is for instance 

\ {gV ~ gV) [n-p]=Au- Ad; \ (gf~ + gf~) [n - p) = Au - Ad; (2.25) 

I (gV~ - gf~) \p] - \ (gV + - gT) [»] = a, - Ac; 

\ (gV' + gf) \P] ~ I {gV + + gf + ) [»] = A," - Ac; (2.26) 
i (gY + - gf + ) \p + n]=Au + Ad + 2As; 

\ (gV~ +gf~)[p + n}= Au + Ad + 2 As. (2.27) 

The structure functions measured from virtual photon scattering eq. (|2.12j ) do not provide 
any further independent linear combination of parton distributions. 

It follows that a complete separation of the four active flavors and antiflavors at a 
fixed scale above charm threshold from inclusive nucleon structure functions only is not 
possible. However, the separation is possible by comparing structure functions above 
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and below threshold. For instance, neglecting Cabibbo-suppressed terms, below charm 
threshold eq. ( |2.24j) determines the down polarized (quark and antiquark) distributions, 
eq. ( |2.26| ) gives the (intrinsic) charm, and then either eq. ( |2.23| ) or eq. ( |2.25| ) give the up 
distribution, the other equation providing a consistency check. Hence, the up, down and 
(intrinsic) charm can always be determined below threshold. In fact, the deviation from 
zero of eq. (|2.26|) provides a simple way of testing for the presence of intrinsic charm or 



anticharm, at least within the limitations of the leading-order approximation. The strange 
distribution can be determined by comparing e.g. eq. ( |2.24j ) below and above threshold. 

Alternatively, above threshold, individual flavor contributions to the structure func- 
tions can be separated by tagging the final-state quark: e.g. the strange contribution can 
be measured by tagging charm in the final state, as is routinely done in unpolarized ex- 
periments ||. A simpler option consists in assuming that the charm distribution vanishes 
below threshold, and is generated dynamically by perturbative evolution above threshold. 
In such mentioned above, the up and down components can be determined from 

eqs. ( |2.23| ),( |2.24| ) below threshold, and the strange distribution can then be determined 
from eqs. ( p.26|) ,( ^27|) above threshold. 

As already mentioned, the virtual photon structure functions eq. ( |2.12[) do not pro- 



vide any further information. However, they do provide independent consistency checks. 
Specifically, as is well known, the isotriplet combination of gj is particularly interesting 
in that its first moment is determined by the Bjorken sum rule. Now, the same isotriplet 
combination of parton distributions is measured by neutral-current and charged-current 
structure functions: 

QgT \p-n] = [gf + - gf~) \p-n} = (Au + Au) - {Ad + Ad) . (2.28) 

Likewise, the strange distribution can be also determined by combining neutral- and 
charged-current structure functions: 

6^7* \P + n}-i ( 9 V + + gf~) W = {gf + - 9?~) [p + n} = (Ac + Ac) - (As + As) . 

(2.29) 

Finally, notice that a determination of the pure singlet quark distribution AS eq.( |2.14[ ) 
can be directly obtained combining charged-current g\ data: 

(gV + + gf~) [p] = (gV + + gV~) [»] = as+. (2.30) 

This is to be contrasted with the analogous measurement obtained from global fits to gj 
data, which can only be done (for first moments) using information from octet /3-decays 
and SU(3) symmetry, or else (for all moments) by using scaling violations. A potentially 
very accurate measurement of this quantity, whose smallness is at the origin of the so-called 
'nucleon spin crisis', is thus possible. 

Of course, all combinations of parton distributions given in this section are only cor- 
rect at leading order, and may thus only be used for rough estimates. Next-to-leading 
order corrections, which are necessary for a more detailed treatment, would modify these 
relations. In practice, an accurate determination of parton distributions can be achieved 
by means of a global fit, which will include all these relations and their next-to-leading 
order modifications. 



8 



3. Structure functions and parton distributions at next-to-leading order 

Beyond leading order, structure functions are no longer linear combinations of quark 
distributions, rather, they are obtained by convoluting quark and gluon distributions with 
perturbative coefficient functions; furthermore, the splitting functions which govern their 
evolution eq. ( |2.16| ) acquire a dependence on Q 2 through a 8 , and in particular no longer 
satisfy the simple relations eqs. ( |2.2q )-( pT2~l~| ). 

3.1. Next-to-leading coefficient functions 

The general expression of the leading-twist polarized charged-current structure func- 
tions beyond leading order is 

gY ± {x,Q 2 ) = i [AC S+ ® AS+t AC NS - ® {Aq»f ~ + Aq» s ~)] + 2[n f /2]AC g ®Ag 
gf ± (x, Q 2 ) = ~\ [AC£_ ® AE" =f AC* NS+ ® {Aq» d s+ + A^ s +)] , i = 4, 5, 

(3-1) 

where rif is the number of active flavors, and [riff 2] is the integer part of nf/2. In this 
equation, Ag is the polarized gluon distribution, the singlet quark distribution AT, ± and 
nonsinglet distributions Aq^ s± were defined in eq. ( p.!4|) and eq. ( |2.15|) respectively, 
and all coefficient functions admit a perturbative expansion of the form AC(x,a s ) = 
AC (0) + §^C (1) (a;, a s ) + . . . , where AC (0) = 5(1 - x) for all the quark coefficient functions 

in eq. ( |3.1| ), and AC^ = for the gluon coefficient function. Notice that the polarized 
gluon distribution contributes to the structure function gi, but decouples from g 4 and g$ 
by charge conjugation. 

To next-to- leading order, in fact, for each of the three structure functions all quark 
coefficient functions are equal in standard factorization schemes such as MS: 



ACe+^Cx) = ACWs- (1) (x) = AC «( 



x 



AC^ {l \x) = AC* NS + (1 \x) = AQW(x), i = 4, 5 
so that eq. (|3.1|) reduces to the simpler expression 

^r ±;NLO (x,Q 2 ) = AC q ® gY ±,h ° + 2[n//2] AC g ® Ag 
^ ± ' NLO (-,Q 2 )=AC^,r ± ' LO , i = 4,5, 



(3.2) 



(3.3) 



where g w - LO are the leading-order expressions of the various structure functions, given 
by eq. ( |2.10| ) above charm threshold and eq. ( |2.11| ) below charm threshold. Note, however, 



that the simple relations eqs. ( |3.2|) are already violated at order a 2 s in the MS scheme |31~H . 



The next-to-leading coefficient functions for virtual photon DIS structure functions 



have been known for some time |32[]. Because the underlying partonic subpro cesses are 
the same, the g± coefficient functions for charged-current and neutral-current scattering 
in fact coincide: the only difference is in the couplings between gauge bosons and quarks, 
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which is reflected in the different leading-order expressions eq. ( 2.10 ) and eq. ( 2.12 ). As 
is by now well known |]J, the first moment of the singlet polarized quark distributions 
has an 0(1) scheme dependence, due to the fact that at leading order the first moment 
Ag(l) of the gluon distribution evolves as l/a 8 , i.e. a s Ag{l) is scale independent up to 
next-to-leading corrections. This is a consequence [|33| of the anomaly which affects the 
singlet axial current. A further consequence of this is that it is possible to choose the 
factorization scheme in such a way that the first moment of the singlet quark distribution 
is scale independent (to all perturbative orders) , even though this is not the case in the MS 
scheme. This choice is particularly convenient in that it allows a meaningful comparison 
with constituent quark model predictions and expectations, which are scale independent. 
Even after fixing by this requirement the scheme ambiguity on the quark first moment, 
there remains a residual ambiguity on the other moments, which can be fixed in a minimal 
way by requiring the scheme-change matrix from the MS scheme to be ^-independent. 
This defines the so-called Adler-Bardeen (AB) scheme [04 1. 



The first moment of any quark distribution in the MS and AB schemes are related by 
A?P(1,Q 2 ) = A^ AB (1) - ^(1,Q 2 ), (3.4) 

where Ag(l, Q 2 ) = L dxAq(x, Q 2 ). Correspondingly, the first moment AC S (1, a s ) of the 
gluon coefficient function vanishes in the MS scheme, while in the AB scheme it is equal 
to AC^ B (l,a s ) = — j^. Notice that, below charm threshold, scheme invariance of the 

physically observable structure function gY eq. (|3J]) upon this scheme change (and more 
generally upon any scheme change which mixes the quark singlet and gluon distributions) 
is ensured by the fact that the mixing with the gluon of As and Ad contributions combine 
to give one effective active flavor. The explicit expressions of the quark and gluon gi 
coefficient functions are given e.g. in Ref. |]34] . 



The next-to-leading order coefficient functions for the structure functions g^ and g§ 



have been computed in Ref. [22]. Both in the MS and AB schemes, they are simply related 



to the gi quark coefficient function by 



x) 



ACf } (x) = AC« (x) + C F x(l + 
ACf } (x) = AC^\x) + C F x(l - x) 



(3.5) 



N 2 -l 



where Cf = on • ^ ne ^ ac ^ that these coefficient functions are unequal to each other 
implies in particular that the Callan-Gross-like relation, eq. ( |2.9|) , is violated beyond lead- 
ing order, as expected since the helicity conservation on which it was based only holds at 
leading order. In fact, the relations between the three quark coefficient functions eq. ( |3.5| ) 
are the same as the corresponding relations between unpolarized quark coefficient func- 
tions which are obtained by the replacements eq. ( |2.7|) , again due to the fact that the 
corresponding tensor structures in eq. ( |2.1| ) are the same. However, by charge conjugation, 
there is no gluon contribution to either g^ or g$, so the violation of eq. (|2.9|) is entirely 
given by eq. ( |3.5| ), contrary to the unpolarized case where the Callan-Gross relation is also 
violated by a gluon contribution. Hence, unlike F± and F2 which, beyond leading order, 
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Figure 1: The lowest-order diagram which contributes to P q q and Pq q D ■ 



measure different combinations of parton distributions (so in particular their comparison 
allows a direct extraction of the gluon distribution), and g§ measure the same combina- 
tion of parton distributions, and thus their separate determination does not provide any 
extra information. Therefore, we will henceforth only discuss the structure function #5, 
since #4 can be entirely determined from it. Once again, this fact can be tested directly on 
the data, by performing a three-parameter (i.e. gi, g±, and #5) fit of the y distributions. 

3.2. Evolution equations 

Beyond leading order the structure of the evolution equations is complicated |30 



by the fact that now both P q q and Pj^q., as well as their polarized counterparts AP g(? 
and AP^® no longer vanish: at next-to-leading order they are all given by the diagram 

displayed in Fig. 1. Furthermore, already at NLO P^ D 7^ P qq , because [^5| there are two 
identical quarks in the final state in the diagonal transition, but not in the non-diagonal 
one. However, 

jyND, NLO tjND, NLO /o c\ 

r qq — *qq 1 K^- ) 

because at this order both processes proceed through the diagram of Fig. 1. 

Coming to the polarized case, it is useful to recall that unpolarized and polarized 
splitting functions are respectively |29| the sum and difference of helicity conserving and 



helicity flipping transitions: 

p = pTT + pti j AP = P TT -P U , (3.7) 

where the arrows here refer to the helicities of the parton before and after the splitting. 
Due to the collinearity of the splitting, parallel (antiparallel) helicities are equivalent to 
parallel (antiparallel) spins. Because the quark-gluon coupling conserves helicity, it follows 
that in the helicity-flipping transition there are two identical particles in the final state 
(including the spin quantum number), but not in the helicity conserving one. Therefore, 

P? 9 U *P? 9 DU , (3-8) 
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but not for the opposite spin configuration: in fact, 

pDTT,NLO _ piVDTT.NLO. 

pJVDTT.NLO pJVDU, NLO pND TT, NLO pJVD TJ.,NLO. /o in\ 

99 99 99 99 ' ^O.lU^/ 

pND t|, NLO _ pD tli NLO 1; Q 

where the second equation follows from the fact that the equality of splitting functions 
eq. (|3.6|) holds for any spin configuration, and the last equation is a consequence of the 
fact that, because of helicity conservation, the flavor-diagonal but spin-flip splitting also 
proceeds only through the diagram displayed in Fig. 1. 

Now, eq. ( j3.8|) immediately implies that the C-even and C-odd nonsinglet splitting 
functions are not the same: AP£ S 7^ AP^ S , like [^5| their unpolarized counterparts. 



However, using eq. ( |3.1U| ) in the explicit expression for the nonsinglet splitting functions 
eq. ( |2.19| ) it immediately follows that 

AP~f h0 = AP~> NLO = AP q D q NLO - AP%> NLO = AP~> NLO (3.12) 

and similarly for their unpolarized counterparts (with AP — > P). Namely, at next-to- 
leading order there is only one C-odd polarized (and one unpolarized) splitting function, 
and thus each C-odd distribution Aq~ eq. (|2.13| ) for the z-th flavor evolves independently: 

d A oi s (t) A „ . 
The difference between singlet and nonsinglet C-odd splitting functions starts at order a s , 



but is expected to be very small |36| at that order. 

Therefore, at next-to-leading order the full set of evolution equations is given by sup- 
plementing the standard singlet and nonsinglet C-even evolution equations eq. (|2.16|) with 



the C-odd evolution equation eq. (|3.13|) . The two- loop polarized splitting functions have 



been computed recently pi|, and are specifically discussed in the AB and related factor- 
ization schemes in Ref. [[34]. The C-odd splitting function, and in fact all the nonsinglet 
splitting functions [j3?| , can be determined from their unpolarized counterparts. Indeed, 
using respectively eq. ( |3.9|) and ( |3.11|) in the definitions of eq. fl2.7Q of AP^ S and AP^ S 
we get 

a p>NS, NLO _ pNSfl, NLO _ pNS, NLO 

99 99 99 ( } 

\pN_S, NLO _piVStl,NLO _ p N_S, NLO K^-^) 

99 _ 99 _ 99 ' 



which imply 
In Mellin space, the anomalous dimensions 



AP^ NLO = P^ NLO (3.15) 



1NS (N)= f dxx N - 1 P NS (x) (3.16) 

do not admit an analytic continuation for all N; however, their even and odd moments can 
be analytically continued separately for all N. Because even (odd) moments correspond 
to C-even (C-odd) operators, the respective analytic continuations give the moments of 
P+ s = AP~ S and P^ s = AP+ S Equations (|3~9|) -(|3~Tl|) and their consequences hold 



in the MS scheme, and remain valid in schemes, such as the AB scheme |34| where the 
nonsinglet and C-odd quark anomalous dimensions are the same as in MS. 
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4. Theoretical expectations and constraints 



Even though neutral-current structure function data do not give any information on 
the relative size of quark and antiquark (or 'valence' and 'sea') distributions, some exper- 



imental information is available from semi-inclusive polarized experiments []38]. However, 
the only conclusion which can be drawn from these experiments at present is that the 
first moments of the Au(x) and Ad(x) distributions are much smaller than the C-even 
combinations Au(x) + Au(x) and Ad(x) + Ad(x), and in fact compatible with zero. When 
the same data are folded in a global NLO analysis |3!| they lead to an indication that the 
first moment of the Au(x) distribution is positive, and roughly 10% of the first moment 
of Au(x) + Au(x), but since this is smaller than the typical systematic uncertainties on 
polarized first moments the only solid conclusion that can be drawn from this is, again, 
that the light antiquark distributions are significantly smaller than the C-even combina- 
tions measured in inclusive experiments. Therefore, a wide variety of theoretical scenarios 
is compatible with present-day data. However, even in the absence of experimental infor- 
mation, some constraints on the relative size of the Aq and Aq distributions are imposed 
by the requirement of positivity of physical cross-sections. 

4-1. Positivity bounds 



Positivity bounds |RJ on polarized structure functions follow from their definition in 



terms of a difference of cross-sections eq. ( |2.6|) , while their unpolarized counterparts are 
defined from the sum of the same cross-sections. Expressing structure functions in terms 
of parton distributions then leads to bounds on the parton distributions themselves. If 
the leading-order expressions of structure functions are used, these bounds coincide with 
those which are obtained in the naive parton model by interpreting parton densities as 
(positive-semidefinite) number densities, while beyond leading order they differ from the 
naive partonic bounds by a calculable amount. Of course, given accurate measurements 
of all the relevant quantities, positivity will take care of itself, i.e. it will just be a trivial 
consequence of the fact that cross-sections are positive. However, in the presence of 
incomplete information, positivity imposes nontrivial bounds on the possible outcomes of 
future experiments. 

In Ref. | 40|1 , positivity bounds on the C-even quark distributions and the gluon were 



derived from their definitions in terms of virtual photon DIS, while it was mentioned 
(but not proven) that charged-current DIS implies similar bounds on quark and antiquark 
distributions separately. In order to prove such bounds, it is enough to consider the cross- 
section asymmetry for gauge boson-hadron scattering Ai, analogous to that discussed in 
Ref. |[40|1 , but with the 7* beam replaced by a VF ± beam. To this purpose, consider the 
cross-section for scattering of a W with helicity Xw off a hadron with helicity A p , given 
by 

a w± (X w ,X p ) =Ke»(\ w y*(\ w )W™ ± (\ p ), (4.1) 

where K contains the coupling constants and the flux and phase-space factors. The polar- 
ization vector is e M (±l) = -^(0, 1, ±i, 0). 



A straightforward calculation using the definition eq. ( |2.1| ) of the hadronic tensor leads 
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to 



a w (X Wl X p )=K (F^-X w F^/2)+^(g^+X w gl 



w- 



(4.2) 



Therefore, we get 



A^(X w ,x,Q' A ) 



a{X w , -1) - a(X w , +1) 
a(X w ,-l) + a(X w ,+l) 
gY ± {x,Q*) + X w gY ± 



x,Q< 



(4.3) 



F™±{x,(p)-\ w F™±{x,(p)/2 



The asymmetry eq. ( |4.3| ) is equal to half (because of the factor discussed after eq. fl2.7p ) the 
ratio of the cross-sections eq. ( |2.6|) and (|2.4| ) evaluated at y = 1. This is a consequence of 
the fact that when y = 1 the lepton and the gauge boson are collinear, and their helicities 
must then be aligned by helicity and momentum conservation, so Xw = A^. 

The leading order positivity bounds now follow immediately by noticing that, because 
of eq. ( |2.10|) and its unpolarized counterpart, the four combinations in the numerator and 
denominator correspond respectively to the sum over generations of the two unpolarized 
and polarized flavors and antiflavors in each generation: e.g., above top threshold 



|< + < LO (- 



-l,x,Q 2 



\Adjx, Q 2 ) + Asjx, Q 2 ) + Abjx, Q< 
d(x,Q 2 ) + s(x,Q 2 ) +b(x,Q 2 ) 



< 1. 



(4.4) 



The three other combinations give the analogous expressions for it-type quarks and for 
the w-type and d-type antiquarks. Because the bound \Ai\ < 1 must be satisfied as a 
matter of principle for any choice of target, below and above each threshold, and also for 
the process where the flavor of the final-state quark is tagged, we immediatley get the 
bounds of Ref. 
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|A^(a;,Q 2 )| < qi (x,Q' 



|A^(x,Q 2 )| < qi(x,Q' 



(4.5) 



for any flavor i. It is easy to see that these conditions are sufficient to ensure positivity of 
the physical lepton-hadron cross-section. Indeed, using the leading-order expression of the 
structure functions eq. (|2.10|) in the cross-section eq. (Ol) > (HI) ' ^ ne positivity conditions 
for the z^-hadron cross section is 



|Au - (y - l) 2 Ad\ <u-{y- l) 2 d, 



(4.6) 



which is manifestly true if eq. ( |4.5|) holds. Similar conditions hold for the other choices of 
lepton beam. 

As discussed in Ref. [|4"0[j , the Altarelli-Parisi equations ( |2.16| ) imply that if the posi- 
tivity bounds are respected at some scale Qq by a given set of polarized and unpolarized 
parton distributions, they will also be respected by the same parton distributions for all 
Q 2 > Qq. Conversely, this implies that the bounds will always be violated at sufficiently 
low scale. However, at low scale the leading-order approximation breaks down, and it 
is necessary to consider higher-order corrections. When the next-to- leading order expres- 
sions of the structure functions eq. (|3.3|) are used, the positivity bounds are best written 
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Figure 2: Leading-order positivity bounds for the C-odd up distribution (solid). The 
Au = curve (dashed) is also shown. 



in terms of the Mellin transforms (defined as in eq. ( |3.16| )) of the coefficient functions and 
parton distributions, so that all convolutions turn into ordinary products. In such case, 
the next-to-leading order bound eq. (|4.4|) , including for simplicity only the lightest flavor 
(as one gets below charm threshold neglecting Cabibbo-suppressed contribution) becomes 



[l + ft (A<#> + ACf } ) 






Au+^AC^Ag 




1 + & + c 











< 1, (4.7) 



where C q , C g and Cs are the unpolarized F\ and F3 coefficient functions, and all parton 
distributions are functions of Q 2 and the Mellin variable N. 

The impact of the gluon correction on the quark positivity bounds has been extensively 



discussed in Ref. [[401] . Because gluons only mix with the C-even combination of quark 
distributions, and thus enter eq. (|4.7| ) through the g\ contribution, the conclusions of that 
reference are unchanged. Namely, because the gluon distribution is peaked at small x, 
the correction to the leading-order positivity bound is only sizable for small moments 
N ~ 1.5, but there the bound itself is very loose because g(N) and q(N) diverge as N — > 1 
while Ag(N) and Aq(N) remain finite, so there is no bound as N — *■ 1. Thus, unless the 
scale is outside the perturbative region (for instance Q 2 1 GeV 2 ), gluon-driven next-to- 
leading order corrections to the positivity bound are essentially negligible. On top of gluon 
corrections, eq. Ql.71) also contains a correction due to the fact that the quark coefficient 
functions for g\ and g$ are not the same, as per eq. ([3.50 . This difference, however, is 
actually quite small and rapidly decreasing at large x: indeed, from eq. ( |3.5| ) we get that 
AC^iN) - AC ( q 1] (N) = ■ Hence even at Q 2 ~ 1 GeV 2 where a s ~ 0.5, and 

even for, say, N = 1.5 the correction is around 1%, and rapidly decreasing as iV increases: 
but, again, at low iV the bound is not relevant because of the vanishing of as iV — > 1. 
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Hence, to very good approximation, the next-to-leading order positivity bounds are 



g(iV) 



- 1 + ^AC( 1 )(A^)' 



i.e. the only significant correction to the leading-order bound is due to the fact that 
the polarized and unpolarized quark coefficient functions are unequal to NLO. Because 

C (1) (A0 

however liniAr^oo A ^ (1)( . jV ^ = 1, this difference is again only significant if N is not too 



large: in practice, the correction was shown in Ref. |1D] to lead to a modification of a few 
percent of the leading-order bound for moments around iV ~ 5, where positivity bounds 
can be relevant. We conclude that for all practical purposes it is enough at the present 
stage to consider leading-order positivity bounds, since next-to-leading corrections, in the 
form of eq. ( f4.8|) , are only relevant at the level of accuracy of a few percent, and then only 
for a limited range of moments, relevant for the description of the shape of the parton 
distributions when x <i 0.1. 

Using the experimentally measured values of the unpolarized parton distributions and 
the polarized C-even distributions, we can then turn the pair of positivity bounds eq. ( f4~5p 
on quark and antiquark distributions for each flavor into a bound on the admissible value 
of the unknown C-odd quark distribution Aq~ eq. ( |2.13| ). The uncertainty on the bound 
will be dominated by, and thus essentially of the same size as, that on the polarized 
quark distributions |5|], which are much more poorly known than the unpolarized ones. 
In Figs. 2-4, we show the positivity bounds on the up, down and strange distributions at 
Q 2 = 5 GeV 2 obtained in this way by using the CTEQ5 |4l]] unpolarized up and down 
quark and antiquark distributions and the BPZ |E2| strange distribution, and the polarized 
parton distributions from Ref. Use of the BPZ fit is motivated by the fact that in this 
fit a particularly accurate determination of the unpolarized strange quark and antiquark 
distributions was achieved by means of a detailed analysis of available neutrino DIS data, 
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Figure 4: Leading-order positivity bounds for the C-odd strange distribution (solid). 
The As = curve (dashed) is also shown. 

whereas a separate determination of the strange quark is not attempted by the CTEQ group 
(which only determines an SU(3) symmetric sea with corrections for the u — d asymmetry). 
In fact, it turns out that the C-even polarized distribution from Ref. [ffl is incompatible 
with the leading-order positivity bound obtained from the CTEQ5 unpolarized strange 
distribution at large x (x ^ 0.5). Even though the violation is within theoretical errors, 
it is suggestive that no such violation is found when the nominally more accurate BPZ 
strange distribution is used. 

In Figs. 2-4, we also display the Aq~ = Aqf , i.e. Aqi = curve for each flavor (dashed 
curves). It is clear that a large value of the antiquark distribution is not compatible with 
positivity for both the up and down quark distributions. Indeed, the curve Aq = Aq (i.e. 
Aq~ = 0) violates the positivity bound for the up quark when x ~ 0.1 and for the down 
quark when x *5 0.2, while when x ~ 0.3 only a very small value of the polarized antiquark 
distribution is allowed (about an order of magnitude smaller than the polarized quark 
distribution). On the other hand, a vanishing light antiquark distribution is compatible 
with positivity for all x. In the case of the strange distribution, instead, both As = and 
As = As (or indeed As = —As) are compatible with positivity. As one should expect, 
for the strange polarized distribution, positivity does not force the antiquark to be smaller 
than the quark distribution. 

4-2. Theoretical scenarios and the spin of the proton 

As we have already mentioned in the introduction, one of the main reasons of interest 
in polarized quark distributions is the unexpected smallness of the nucleon axial charge 
which has been determined in the first generation of polarized DIS experiments. A clarifi- 
cation of the physics behind this requires a determination of the detailed polarized parton 
content of the nucleon. A review of the various explanations which have been proposed for 
this experimental fact is beyond the scope of this work. However, it is useful to sketch var- 
ious scenarios for the polarized content of the nucleon which are representative of possible 
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theoretical alternatives, and which could be tested in future experiments. 



Firstly, it should be noticed that even though current data give a value of the 
axial charge which is compatible with zero, they cannot exclude a value as large as 
ao(10 GeV 2 ) = 0.3 0. Also, the current value is obtained by using information from 
hyperon /3-decays and SU(3) symmetry. Clearly, the theoretical implications of an exact 
zero are quite different from those of a value which is just smaller than expected in quark 
models. It is thus important to have a direct determination of the axial charge. If a small 
value is confirmed, it could be understood as the consequence of a cancellation between 
a large value of the scale-independent (i.e. AB-scheme) quark first moment, and a large 
gluon first moment. Indeed, since in the MS scheme the axial charge and singlet quark 
first moment coincide [|33]], eq. ( |3.4| ) is immediately seen to imply that a large gluon con- 



tribution can lead to a small value of the axial charge even when the AB-scheme quark 
is large. In this ('anomaly') scenario the up, down and strange polarized distributions 
in the AB-scheme are close to their expected quark model values, so in particular the 
strange distribution is much smaller than the up and down distributions. In Ref. |[43|| , 



this cancellation of quark and gluon components has been derived from the topological 
properties of the QCD vacuum ||3| (and thus further predicted to be a universal property 
of all hadrons). 

If instead the polarized gluon distribution is small, the smallness of the singlet axial 
charge can only be explained with a large and negative strange distribution. In this 
case, the scale-independent first moment of the singlet quark distribution is also small. 
This scale-independent suppression of the axial charge might be explained by invoking 



non-perturbative mechanisms based on instanton-like vacuum configuration [J3|. In this 
'instanton' scenario the strange polarized distribution is large and equal to the anti-strange 
distribution, since gluon-induced contributions must come in quark-antiquark pairs. 

Another scenario is possible, where the smallness of the singlet axial charge is due to 
intrinsic strangeness, i.e. the C-even strange combination is large, but the sizes of As and 
As differ significantly from each other. Specifically, it has been suggested that while the 
strange distribution (and specifically its first moment) is large, the antistrange distribution 
is much smaller, and does not significantly contribute to the nucleon axial charge J55J . This 



way of understanding the nucleon spin structure is compatible with Skyrme models of the 
nucleon, and thus we will refer to this as a 'skyrmion' scenario ||. 

Therefore, the main qualitative issues which are relevant for the nucleon spin structure 
are to assess how small the axial charge is, to determine whether the polarized gluon 
distribution is large, and then whether the strange polarized distribution is large, and 
whether the strange polarized quark and antiquark distributions are equal to each other 
or not. More detailed scenarios might then be considered, once the individual quark and 
antiquark distributions have all been accurately determined. For instance, while the up 
and down antiquark distributions are small, they need not be zero, and in fact they could 



be different from each other [^6|, just like their unpolarized counterparts appear to be. 
Investigating these issues could shed further light on the detailed structure of polarized 
nucleons. 
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Figure 5: Charged-current event rates, in units of 10 6 , as a function of Lab- frame 
neutrino spectra, for several detector and beam configurations. The dashed lines on the 
left include cuts on the final-state muon [E^ > 3 GeV) and on the final-state hadronic 
energy {Ehad > 1 GeV). The solid lines have no energy-threshold cuts applied. The 
three set of curves correspond to different detector radii (50, 10 and 5 cm, from top to 
bottom). 

5. Polarized neutrino-DIS at the front-end of a muon ring 

The parameters of a realistic neutrino-factory complex are still under active study, and 
it is expected that the accelerator configuration will evolve and upgrade in both energy and 
luminosity after the beginning of operations |37]]. For the present study we shall assume 
the following beam parameters, which are considered as typical of the initial operations 
of a neutrino factory: muon beam energy, = 50 GeV; length of the straight section, 
L = 100 meters; distance of the detector from the end of the straight section, d = 30 meters; 
number of muon decays per year along the straight section, = 10 20 ; muon beam angular 
divergence, 0.1 xm^/E^; muon beam transverse size a x = a y = 1.2 mm. While the concept 
of "conservative" cannot be applied to any of the projections for the neutrino factory, it is 
nevertheless reasonable to expect that should one ever be built, its performance would not 
be inferior to what assumed here. We also assume a cylindrical detector, with azimuthal 
symmetry around the beam axis, parametrized by its radius R = 50 cm and a thickness 
of 10 gr/cm 2 . The collected statistics scales linearly with the detector thickness, while the 
dependence on other parameters such as the radius or the length of the straight section is 
clearly more complex. 

The neutrino spectra are calculated using standard expressions for the muon decays 
(see e.g. Ref. |15|]). For simplicity we shall confine ourselves to the case of and charged 
current DIS. The laboratory-frame neutrino spectra, convoluted with the charged-current 
interaction cross-sections, are shown for several detector and beam configurations in Fig. 5 
(Efj, = 50 GeV) and in Fig. 6 (E^ = 100 GeV). 
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Figure 6: Same as Fig. 5, but for 100 GeV muon beams. 



To determine the statistical accuracy with which the individual structure functions 
and their flavor components are expected to be measured we shall exploit the different 
y dependence of the separate F{ and components of the cross-section. An important 
feature of the neutrino beams from muon decays is their wide-band nature. This allows to 
modulate the y dependence for fixed values of x and Q 2 using the neutrino energy: 

O 2 

y = -^r- (5.1) 

y 2xmE v y ' 

The separate measurement of the recoil muon energy and direction, and of the recoil 
hadronic energy, enable the event-by-event extraction of x,y and Q 2 . We assume perfect 
experimental resolution on the determination of these quantities, after imposing a cut on 
the minimum energy of the muon (> 3 GeV) and of the hadronic system (> 1 GeV). 
Since we are just interested in predicting the statistical accuracy in the determination of 
the structure functions, we also assume for simplicity the absence of power-suppressed 
corrections, and impose the Callan-Gross relation. Given the large statistics expected at 
the neutrino factory, it will nevertheless be possible to use the data to extract the higher- 
twist components and to separate the individual contributions of F\ and F 2 (#4 and #5), 
without any theoretical assumption. This issue is studied in some detail for the case of 
unpolarized distributions in Ref. ]16[ . 



We thus produced y distributions by generating events within different bins of x and 
Q 2 , and performed minimum-x 2 fits of the generated data using the cross-section eq. ( |2.4f) . 
For each bin, the values of x and Q 2 at which we quote the results are obtained from the 
weighted average of the event rate. As an input, we used the CTEQ4D set of parton 



distributions [48 . The dependence on the parameterization of the parton distributions 



is very small, and will be neglected here. We verified that other recent sets of parton 
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Figure 7: x and Q 2 binning for the generation of charged-current events. The crosses 
correspond to the weighted bin centers. 



distributions give similar results. The fit to the y distributions at fixed x and Q 2 for a 
fully polarized target gives then the value of the combinations F<i ± 2x^5 and F% ± 2g\. 

The absolute number of events expected in each bin is scaled by the total number 
of muon decays; this number of events determines the statistical error on the individual 
structure functions obtained through the fit. Polarization asymmetries are extracted by 
combining data sets obtained using targets with different orientations of the polarization. 
The statistical accuracies with which the combinations can be performed depend on the 
statistical content of each individual data set. Since the polarization asymmetries are 
small relative to the unpolarized cross-sections, the absolute statistical uncertainties on 
the extraction of polarized structure functions will have a very mild dependence on the 
value of the polarized structure functions themselves; they will be mostly determined by 
the value of the unpolarized structure functions (which to first approximation fix the overall 
event rate), and by the polarization properties of the target. 

For simplicity, we therefore calculate directly the expected statistical errors o~f 2 ,f 3 on 
the extraction of F2 and F3 using unpolarized targets, and relate them to the errors on 



the polarized cross sections using the following relation given in [17 



a gi =F%V2 aiJ ^, (5.2) 

where aij = 1 for = (1,3) and o^j = 1/x for = (5,2), and where F v % is 

a correction factor (always larger than one) which accounts for the ratio of the target 
densities to H2 or D2, for the incomplete target polarization, and for the dilution factor 
of the target, namely the v (or v) cross-section weighted ratio of the polarized nucleon to 
total nucleon content of the target. The factor of y/2 in the numerator reflects the need to 
subtract the measurements with opposite target polarization. 

We generate events in the (x, Q 2 ) bins shown in Fig. 7. Twenty equally-spaced bins 
in the < y < 1 range are used for the y fit. The total number of x bins varies in 
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Figure 8: Absolute errors on proton structure functions for the bins of Fig. 7. 



different Q 2 bins because of kinematic acceptance and minimum energy cuts. The values 
of the uncertainties in the determination of the eight charged-current structure functions 
(<7i and g$ with the two available beams and targets) are assigned at the cross-section 
weighted bin centers. We then get the absolute errors on the structure functions displayed 
in Figs. 8-9 for proton and deuterium targets respectively, using the p-butanol and D- 
butanol target [f|9| correction factors given in Ref. [fl7| , namely = 2.6, = 1.6 and 
Fj? D = 4.4 (for a more complete discussion of polarized targets and their complementary 
properties, see [17] and references therein). We have assumed a luminosity of 10 20 muons 



decaying in the straight section of the muon ring for each charge, for each target, and for 
each polarization. Assuming that only one polarization and one target can run at the same 
time, this means eight years of run. While the number of muons may not be dramatically 
increased, the integration time can be reduced by a large factor if the target thickness can 
be increased over the conservative 10 gr/cm 2 assumed here, or if different targets can be 
run simultaneously. 



6. Determination of polarized structure functions at a neutrino factory 

We can now study how charged-current DIS data may be used to determine the 
polarized parton content of the nucleon. The purpose of this exercise is twofold. First, this 
will give us a first indication of the expected shape and scale-dependence of the charged- 
current structure functions. Furthermore, on a more quantitative level, we will be able 
to assess the level of accuracy with which charged-current DIS experiments will allow a 
determination of individual flavors and antiflavors. Of course, both these tasks can only 
be achieved by making assumptions on the expected flavor content of the nucleon, and on 
the expected experimental accuracy. To this purpose, we consider the theoretical scenarios 
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Figure 9: Absolute errors on deuteron structure functions for the bins of Fig. 7. 



which we have discussed in Sect. 4, and we assume the availability of data such as discussed 
in Sect. 5. 

In order to implement the theoretical scenarios of Sect. 4, we need to take as a starting 
point a set of parton distributions which summarize the current knowledge of the polarized 
structure of the nucleon. To this purpose, we adopt the type-A fit of Ref. f§], which is 
defined as follows. The C-even polarized quark distributions AE + and 



Asa = A^/+ 
Aq 8 = Aq» S + + Aq^ + 



(6.1) 



ds 

2 _ i n„\r2 



and the polarized gluon distribution Ag at the initial scale Qq = 1 GeV are all taken to 
be of the form 

Af{x,Ql)=Ar f 7i f x a t(l-xf'{l + >y f x 5 '), (6.2) 

where the factor A// is such that the parameter r/f is the first moment of Af at the initial 
scale. The nonsinglet quark distributions A^3 and Ag§ are assumed to have the same x 
dependence, while the parameter rjs, corresponding to the first moment of Aqg, is fixed to 
the value r]s = 0.579 from octet baryon decay rates using SU(3) symmetry. Furthermore, 
7e = 7 g = 10, ds = ds = 0.75, Sj; = S g = 1. All other parameters in eq. ( |6.2|) are 
determined by the fitting procedure. 

This set of polarized parton distributions is particularly useful for our present purposes 
because a detailed analysis of statistical and systematic uncertainties on the first moments 
based on it was performed in Ref. Q . The set A is adopted because all parton distributions 
in this set are amply within positivity bounds [40|. However, some new data [^0| from the 



E155 collaboration, and the final data set from the SMC collaboration have appeared since 
the publication of Ref. H . Therefore, we have repeated the fit including these new data, 
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in order to ensure that no available information is neglected. This gives us a data set of 
176 neutral-current data points. 

The best-fit values of the first moments of the C-even parton distributions at the 
initial low scale Qq = 1 GeV 2 thus obtained are listed in the first column of Table 1, 
together with the errors from the fitting procedure. The fit is performed, and all results 
are given, in the AB scheme. These values are not far from those of Ref. M. They are 
further affected by theoretical uncertainties, for which we refer to The detailed shapes 
of parton distributions are similar to those of Ref. 0], and likewise not very accurately 
known. As mentioned above, the value of r]g in this fit is fixed using the SU(3) flavor 
symmetry and the measured octet baryon decay constants. Even though the nominal 
uncertainty on this SU(3) value of rjs is relatively small (around 5%), the uncertainty due 
to failure of exact SU(3) symmetry could be as large as 30% [|51|. This effect is taken 
into account in the estimated theoretical uncertainties given in Ref. |4]]. The subsequent 
three rows of the table give the values of the first moments of Aq(x,Q 2 ) + Aq(x,Q 2 ) at 
Q 2 = 1 GeV 2 for up, down, and strange, obtained combining the singlet and nonsinglet 
quark first moments above. Finally, we give in the last row the value of the singlet axial 
charge ao at the scale Q 2 = 10 GeV 2 . 

Because the first moment of the gluon distribution in this fit is quite large, we can 
take this global fit as representative of the 'anomaly' scenario discussed in Sect. 4.2, even 
though the strange distribution is not quite zero. In order to construct parton distributions 
corresponding to the other scenarios discussed in Sect. 4.2, we have also repeated this fit 
with the gluon distribution forced to vanish at the initial scale. This possibility is in fact 
disfavored by several standard deviations; however, once theoretical uncertainties are taken 
into account a vanishing gluon distribution can only be excluded at about two standard 
deviations [|J, and thus this possibility cannot be ruled out on the basis of present data. 
The results of this fit for the various first moments are displayed in the second column of 
Table 1. 

We can now use these parton distributions to construct the unknown C-odd parton 
distributions. We construct three sets of parton distributions, corresponding to the three 
scenarios of Sect. 4.2. In all cases, we assume Au(x) = Ad(x) = 0. Furthermore, as 
the 'anomaly' set we take the 'generic' fit of Table 1 with the assumption As(x) = 0, 
the strange distribution for this set being relatively small anyway. As 'instanton' and 
'skyrmion' parton sets we take the Ag = fit of Table 1, with As = As in the former case, 
and As = in the latter case. With these choices all quark and antiquark distributions 
are fixed, and thus all structure functions can be computed. 

We generate for each of these three scenarios a set of pseudo-data, by assuming the 
availability of neutrino and antineutrino beams, and proton and deuteron targets, in the 
(x, Q 2 ) bins of Fig. 7 and with the errors displayed in Figs. 8,9. Although experimentally 
event rates remain sizable even in the large-x region, we do not include data with x > 0.7 
because, in this region, leading twist next-to-leading order perturbation theory is not 
reliable. We discard data points whose uncertainty is larger than 50, the typical size of 
the structure functions being of order one. We finally use these uncertainties to generate 
data gaussianly distributed about the values of the structure functions at each data point 
in the three scenarios. We obtain in this way approximately 70 data points for each of the 
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par. 


generic fit 


Ag = fit 


'anomaly' refit 


'instanton' refit 


'skyrmion' refit 


1 ^ 


0.38 ±0.03 


0.31 ±0.01 


0.39 ±0.01 


0.321 ±0.006 


0.324 ± 0.008 


Tin 

19 


0.79 ±0.19 





0.86 ±0.10 


0.20 ±0.06 


0.24 ±0.08 


V3 


1.110 ±0.043 


1.039 ±0.029 


1.097 ±0.006 


1.052 ±0.013 


1.066 ±0.014 


V8 


0.579 


0.579 


0.557 ±0.011 


0.572 ±0.013 


0.580 ±0.012 


Vu 


0.777 


0.719 


0.764 ±0.006 


0.722 ±0.010 


0.728 ± 0.009 


Vd 


-0.333 


-0.321 


-0.320 ±0.008 


-0.320 ±0.009 


-0.325 ± 0.009 




-0.067 


-0.090 


-0.075 ± 0.008 


-0.007 ±0.007 


-0.106 ±0.008 


a 


0.183 ±0.030 


0.284 ±0.012 


0.183 ±0.013 


0.255 ±0.006 


0.250 ±0.007 



Table 1: Best-fit values of the first moments for the data and pseudodata fits discussed in 
the text. 

eight charged-current structure functions. 

We proceed to fit a global set of data which includes the original neutral-current 
data as well as the generated charged-current data. We assign to the generated data the 
estimated statistical errors, and fit including statistical errors only. We do not attempt an 
estimate of the experimental systematic uncertainties, which are very difficult to anticipate: 
our results are meant to provide a benchmark for the best possible situation of negligible 
systematics, and set a target of performance for the planned detectors. The errors assigned 
to the neutral-current data are instead obtained as in our original fits, by adding in 
quadrature the statistical and systematic errors given by the various experimental groups: 
a fit including statistical errors only to these data would be very hard to achieve, and 
somewhat misleading, since the data themselves are affected by the systematics, which 
cannot be assumed to be absent. Since the statistics is dominated by the generated data, 
the final errors on the various parameters should be taken as estimates of statistical errors 
only. 

The fits are performed adopting the same functional form and parameters as in the 
original fit for the C-even parton distributions, except that the normalization of the octet 
C-even distribution ryg is now also fitted. For the C-odd parton distributions, we add six 
new parameters, namely the normalizations of the up, down and strange C-odd distribu- 
tions, and three small-x exponents a (corresponding to an x a small-x behaviour). The 
shape is otherwise taken to be the same as that of the C-even quark distributions. The 
exponents a are included among the refit parameters in order for the refitting procedure 
to be more realistic, in view of the fact that the generated data are distributed around 
the best fit with appropriate errors, but do not coincide with it. The charm distribution 
is assumed to vanish below threshold, and to be generated dynamically by perturbative 
evolution above threshold. 

The best-fit values of all the normalization parameters are shown in the last three 
columns of Table 1, where the rows labelled rj u , rid and r] s now give the best-fit values and 
errors on the first moments of Aq~ for up, down and strange. Comparison of these values 
with those of our original fits leads to an assessment of the impact of charged-current data 
on our knowledge of the polarized parton content of the nucleon. 
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First, we see that the improvement in the determination of the polarized gluon distri- 
bution is minor. This is due to the fact that the gluon distribution is determined by scaling 
violations, and thus a precise determination is only possible using data which cover a wide 
range of values of Q 2 . Lacking this, the availability of charged-current data per se does not 
help. Furthermore, since the gluon decouples from g$, the corresponding data do not have 
any effect on the determination of the gluon distribution. In the 'skyrmion' and 'instanton' 
fits, where the charged-current data have been generated assuming a vanishing polarized 
gluon component, the refitted gluon is nevertheless nonzero, though significantly smaller 
than in the 'anomaly' fit. This is due to the fact that the currently available (neutral- 
current) data are also included in the refit, and underlines the fact that these data support 
a gluon distribution which differs significantly from zero. The charged-current data can 
only have a limited impact in modifying this conclusion for the same reason why they lead 
to a modest improvement in the precision of the determination of the gluon. 

Let us now consider the C-even quark distributions. It is immediately clear that the 
precision on the singlet quark first moment is very significantly improved by the charged- 
current data: the error on the first moment of AE + is now of a few percent in comparison 
to about 10% with neutral-current DIS. This follows from the fact that, up to subhead- 
ing corrections, the singlet quark is directly measured by the combination eq. (|2.30| ) of 
charged-current structure functions. This improvement is especially significant since the 
determination of 77s no longer requires knowledge of the SU(3) octet component, unlike 
that from neutral-current DIS, and it is thus not affected by the corresponding theoretical 
uncertainty. With this accuracy, the relatively larger value of the singlet quark component 
found in the 'anomaly' scenario (in the AB scheme) can be experimentally distinguished 
from the smaller value found in other scenarios, at the level of several standard deviations. 
In other words, thanks to the charged-current data, it is possible to experimentally refute 
or confirm the anomaly scenario by testing the size of the scale-independent singlet quark 
first moment. Correspondingly, the improvement in knowledge of the gluon first moment, 
although modest, is sufficient to distinguish between the two scenarios. 

The determination of the singlet axial charge is improved by an amount comparable 
to the improvement in the determination of the singlet quark first moment. Its vanishing 
could thus be established at the level of a few percent. The determination of the isotriplet 
axial charge is also significantly improved: the improvement is comparable to that on the 
singlet quark, and due to the availability of the triplet combination of charged-current 
structure functions eq. (|2.28[) . This would allow an extremely precise test of the Bjorken 
sum rule, and accordingly a very precise determination of the strong coupling. Finally, 
the octet C-even component is now also determined with an uncertainty of a few percent. 
Therefore, the strange C-even component can be determined with an accuracy which is 
better than 10%. Comparing this direct determination of the octet axial charge to the 
value obtained from baryon decays would allow a test of different existing models of SU(3) 



violation [51 



Coming now to the hitherto unknown C-odd quark distributions, we see that the up 
and down C-odd components can be determined at the level of few percent. This accuracy 
is just sufficient to establish whether the up and down antiquark distributions, which are 
constrained by positivity to be quite small, differ from zero, and whether they are equal 
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Figure 10: The structure functions gY* , (top) and gY , g§ (bottom) for a 

proton target. 



to each other or not. Furthermore, the strange C-odd component can be determined 
at a level of about 10%, sufficient to test for intrinsic strangeness, i.e. whether the C- 
odd component is closer in size to zero or to the C-even component. The 'instanton' 
and 'skyrmion' scenarios can thus also be distinguished at the level of several standard 
deviations. 

Of course, only experimental errors have been considered so far. In Ref. it has 
been shown that theoretical uncertainties on first moments are dominated by the small-a; 
extrapolation and higher-order corrections. The error due to the small-x extrapolation 
is a consequence of the limited kinematic coverage at small x. This will only be reduced 
once beam energies higher than envisaged in this paper will be achieved; otherwise, this 
uncertainty could become the dominant one and hamper an accurate determination of 
first moments. On the other hand, the error due to higher order corrections could be 
reduced, since it is essentially related to the fact that available neutral-current data must 



27 



0.2 



0.1 



0.0 7='. 



x -0.1 



- dotdashed: Q =1 GeV* 

- dotted: Q 2 = 20 GeV 2 



-0.2 



°: gi(x,Q 2 ) 
X: g 5 (x,Q 2 ) 




0.05 



0.1 



0.5 



0.2 



0.1 



x -0.1 



-0.2 



1 




O: 
x : 


g,(x,Q 2 ) 
g 5 (x,Q 2 ) 








-j- ' ' % 








\ \ 
■m ■ 

\ - 
/ - 








/ 

y / 


_ solid: Q 2 = Q 2 (data) 

- dotdashed: Q 2 =l GeV 2 

- dotted: Q 2 = 20 GeV 2 

1 









0.05 



0.1 



0.5 



Figure 11: The structure functions gY + , gV (top) and gY , (bottom) for a 

deuterium target. 



be evolved to a common scale, and also errors are amplified when extracting the singlet 
component from neutral-current data because of the need to take linear combinations of 
structure functions. Neither of these procedures is necessary if charged-current data with 
the kinematic coverage considered here are available. 

The best-fit structure functions corresponding to the 'anomaly' refit (third column of 
Table 1) are displayed as functions of x at the scale corresponding to the bin 4 GeV 2 < 
Q 2 < 8 GeV 2 , and compared with the data in Figs. 10-11 (a few data points with large 
error bars, although included in the fits, are not shown in the figures). Notice that the 
structure functions g\ and g*, always have opposite signs because the (dominant) quark 
component in g\ and g§ has the opposite sign, while the antiquark component has the 
same sign. For comparison, we also display the structure functions at the initial scale of 
the fits, and at a high scale. The good quality of the fits is apparent from these plots. 

Given the poor quality of current knowledge of the shape of polarized parton distri- 
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Figure 12: The combinations of structure functions of eq. ( |2,23| ), and the corresponding 
parton distributions. 



butions, it is difficult to envisage detailed scenarios and perform a quantitative analysis 
of the various shape parameter, as we did for first moments. However, it is possible to 
get a rough estimate of the impact of charged-current data on our knowledge of the x 
dependence of individual parton distributions by considering the combinations of struc- 
ture functions given in eqs. (|2.23|) ,( ^.24f) , which, at leading order, are directly related to 
individual parton distributions. In Figs. 12,13 we show, respectively, the combinations 
of eq. (|2.23j ) and ( |2.24|) for a proton target, together with the pseudodata for the same 
combinations of structure functions, in the bin 4 GeV 2 < Q 2 < 8 GeV 2 . In each figure we 
also display the two parton distributions which contribute at leading order to the relevant 
combination of structure functions at Q 2 = 7 GeV 2 , as well as (a s /7r)Ag at the same scale. 

Let us consider the upper plot in Fig. 12. It is apparent that the expected statistical 
accuracy is very good for all data with x > 0.1. This suggests that an accurate determina- 
tion of the shape of Au+ Ac is possible. Furthermore, it is also clear that Ac (dotted curve) 
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is extremely small compared to Au (solid curve). However, we observe that the difference 
between the Au distribution (solid) and the data is of the order of 15% to 20% for all x 
below 0.4. This difference is entirely due to next-to-leading corrections. Specifically, the 
gluon contribution (dot-dashed curve), which in the AB scheme spoils the leading order 
identification of the quark parton distribution with the structure function [see eq. (|3~i])1 , is 
small but non negligible. Because the various contributions to next-to-leading corrections 
(in particular the gluon distribution) are affected by sizable theoretical uncertainties f|], 
this implies that Au can only be determined with an error which is considerably larger 
than the experimental one. At larger scales, one expects the subleading corrections to co- 
efficient functions to be smaller and smaller, while a residual gluon contribution persists, 
because of the axial anomaly []33| . 

A similar analysis of the lower plot of Fig. 12 tells us that a determination of the shape 
of Au is essentially impossible. This combination of structure functions is the preferred 
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one for a determination of the charm distribution, since perturbatively we expect Ac = Ac, 
and Au is much smaller than Au. Nevertheless, it is apparent from this figure that even 
in this determination of the charm distribution is out of reach. 

A study of the down quark and antiquark distributions can be similarly performed 
looking at Fig. 13. The conclusion for Ad is similar, although perhaps slightly less op- 
timistic, to that for Au: a reasonable determination of its shape is possible, but with 
sizable theoretical uncertainties. The lower plot shows that no significant information on 
the shape of Ad can be obtained from this analysis. 

7. Conclusions and Outlook 

We presented in this paper a self-contained review of the next-to-leading order for- 
malism for the description of charged-current polarized DIS, discussing in particular the 
impact of next-to-leading order corrections to the evolution equations and to the relations 
between polarized structure functions and polarized parton densities. We discussed the 
theoretical constraints on individual quark and antiquark polarized distributions emerging 
from positivity requirements of physical cross-sections. We found that Au(x) and Ad(x) 
are constrained to be much smaller than the valence polarizations Au(x) and Ad(x), re- 
spectively; As(x) and As(x) are instead allowed to have similar size. 

We then focused the phenomenological applications to the case of neutrino scattering, 
which has recently become an intriguing possibility in the context of the neutrino factory 
facilities being explored world-wide. After evaluating the statistical uncertainties which 
are expected to be achievable in the measurement of polarized structure functions, we 
studied the impact of such accuracies on the theoretical models describing the structure 
of the proton spin. In the case of the C-even distributions, our results indicate that the 
singlet, triplet and octet axial charges can be measured with accuracies which are up to one 
order of magnitude better than the current uncertainties. In particular, the improvement 
in the determination of the singlet axial charge would allow a definitive confirmation or 
refutation of the anomaly scenario compared to the 'instanton' or 'skyrmion' scenarios, at 
least if the theoretical uncertainty originating from the small-cc extrapolation can be kept 
under control. The measurement of the octet axial charge with a few percent uncertainty 
will allow a determination of the strange contribution to the proton spin better than 
10%, and allow stringent tests of models of SU(3) violation when compared to the direct 
determination from hyperon decays. 

In the case of C-odd distributions, the up and down components can be determined at 
the level of few percent, allowing the measurement of an antiup and antidown polarization 
at levels small enough to be still compatible with the positivity constraints. The strange 
C-odd component can be measured at the level of 10%, sufficient to test for instrinsic 
strangeness, and thus distinguish between 'skyrmion' and 'instanton' scenarios at a level 
of several standard deviations. 

We also studied the prospects for the determination of the shape of the polarized 
distributions. If one were to assume the leading-order relation between structure functions 
and polarized parton densities, the statistical power of the neutrino factory data would pin 
down the shape of Au(x) with a few-percent precision. Most of this accuracy is however 
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lost when next-to-leading order corrections are taken into account, which mix the quark 
and gluon distributions. The impact of the gluon distribution on the extraction of the 
Am(x) shape is at the level of 15-20%, and any uncertainty on Ag(x) will be reflected on 
this extraction accordingly. The results for Ad(x) are similar to those for Am(j), while 
no significant shape information can be obtained for the sea distributions. In the case of 
As(x) one may be able to use semi-inclusive measurements, with tagged charm quarks in 
the final state, but this study remains to to be done. 

Our analysis did not try to incorporate any estimate of the systematic experimental 
uncertainties, nor of the theoretical systematics induced by the separation of higher-order 
and higher-twist contributions, and small-x extrapolation. Likewise, we did not try to 
include the new information on the spin structure of the proton and on polarized parton 
densities which might become available through future experiments at CERN, DESY and 
RHIC. Our study confirms nevertheless the expectation that polarized DIS experiments 
at a neutrino factory will provide invaluable information on the structure of the proton, 
provided the experimental systematics will be able to match the statistical one. At the 
same time, our work stresses the known fact that QCD corrections to the parton-level 
picture are large, and that a neutrino factory alone may not be able to fully disentangle 
the shape of individual quark flavors unless a firmer knowledge of the polarized gluon is 
achieved. 
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